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Abstract. The possibility that population III stars have reionized the Universe at redshifts greater than 6 has recently gained 
momentum with WMAP polarization results. Here we analyse the role of early dust produced by these stars and ejected into 
the intergalactic medium. We show that this dust, heated by the radiation from the same population III stars, produces a submil- 
limetre excess. The electromagnetic spectrum of this excess could account for a significant fraction of the FIRAS (Far Infrared 
Absolute Spectrophotometer) cosmic far infrared background above 700 micron. This spectrum, a primary anisotropy (AT) 
spectrum times the v 2 dust emissivity law, peaking in the submillimetre domain around 750 micron, is generic and does not 
depend on other detailed dust properties. Arcminute-scale anisotropies, coming from inhomogeneities in this early dust, could 
be detected by future submillimetre experiments such as Planck HFI. 
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1. Introduction 

More accurate measurements of the cosmic microwave back- 
ground (CMB) implies a need for a better understanding of the 
different foregrounds. We study the impact of dust in the very 
earl y universe 5 < z < 15. WMAP data on the CMB polariza- 
tion, iKo gut et alJ j2003) provides a strong evidence for a rather 
large Thomson opacity during the reionization of the Universe: 
r e = 0.17 ± 0.04 (68%C.L.). Although the me chanism of pro- 
ducin g such an opacity is not fully understood. ICen. r] |2002, 
2003) has shown that early, massive population-Ill (Pop III) 
stars could ionize the Universe within 5 < z < 15 (see Fig. [2 
and Fig. [5). Adopting this hypothesis, we discuss the role and 
observability of the dust that is produced by the Pop III stars. 
As we can only conjecture about the physical properties and 
the abundance of this early dust, we adopt a simple dust grain 
model with parameters deduced from the Milky Way situation. 
The dust production is simply linked to the ionizing photon 
production by the stars through their thermal nuclear reactions. 
The low potential well of the small pre-galactic halos allows the 
ejected dust to be widely spread in the intergalactic medium. 
The ionizing and visible photons from the same Pop III stars 
heat this dust. There are no direct measurements of this dust, 
but by means of other results the amount of dust can be es- 
timated. A similar study has been done for a later epoch of 
the un iverse, in which data are more readily available. IPei et al] 
(1999). We use a cosmology with £l tot - Q m + Oa = 1, where 
Q m = Q fc + Q. DM = 0. 1 33 /ft 2 , Cl h = 0.0226 1 h 2 and h = 0.72 as 
advocated by WMAP, Snerg el et all (12003). using WMAP data 
in combination with large scale structure observations (2dF- 
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Fig. 1. Total number of ionizing photons produced from Pop III 
stars per baryon, cf. llCen. Rl2002[ figure 14). The dotted line 
represents a simplified model with a constant photon produc- 
tion, from z = 16, of 8 per unit z per baryon. The results are 
similar. 



Send offprint requests to: Erik Elfgren, e-mail: elf@ludd.luth.se 



GRS + Lyman a). Furthermore, since z » 1 the universe is 
matter-dominated. We relate all cosmological parameters to 
their measurement today so that they have their present-day 
values throughout our calculations. 

We now proceed to compute the abundance and the tem- 
perature of this dust. Consequences on the CMB distortions are 
then to be discussed. 
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Fig. 2. Production rate of ionizing photons from Pop III stars 
per baryon, ^r/n;,. The odd form between each integer z is not 
physical but are due to the fact that the redshift z is a nonlinear 
function of time. 

2. Dust Model 

Here we assume the dust properties to be similar to what is 
observed in our galaxy. For simplicity, we suppose spherical 
dust grains with radius a = 0.1pm and density p g = 2.3 • 10 3 
kg/m 3 .The absorption cross section, cr y , between photons and 
dust can be written as 

o- v = Q v 7ra 2 , (1) 
where we parametrize the frequency dependency as 



q _ | Qof- {jff submm and infra red (IR), 
I 1 visible and ultra violet (UV), 



(2) 



where v r , a r and Qo are normalization constants. There is only 
one indep endent constant wh ich means that we can fix a r = 
0.1 pm. In llDesert et alll990l figure 3) the poorly known knee 
wavelength, A r = c/v r was set to 100 pm. Here, we choose 
40 pm for simplicity, so that early dust radiates mostly in the 
v 2 emissivity regime. Above the characteristic frequency v r the 
spectral index B = \, below (3-2. The exact position of v r 
is not very important for our study because it is mainly above 
the interesting wave-length region ~ 0.3 - 3 mm and it will not 
change the magnitude of the signal. 

In the submm and far IR range, the spectral index is con- 
stant, and with Qo = 0.0088 the assumed opacity agrees well 
with measurements by FIRAS on cirrus cloud s in our galaxy, 
cf.lBoulanger et alJ lll996h : lFixsen et alJ lll998h : lLagache et alJ 
( 1999). In the visible and UV region, the cross section is inde- 
pendent of the frequency because A < 2na. In the submm re- 
gion, the cross section is proportional to the mass of the grain. 

In order to evaluate the significance of the dust during the 
reionization, we calculate the amount of dust present in the 
universe at a given time. The co-moving relative dust density 
is D. d ,o = Pd(z)/((1 + z) 3 p c ), where pd(z) is the dust density, 
z is the red-shift, p c = ^ is the critical density (Hq and G 
are Hubble's and Newton's constants, respectively). The co- 
moving relative dust density as measured today evolves as: 

—^-=J + -J-, (3) 

dz 



where J+ and J are the production and the destruction rate 
respectively. 

The Pop III stars produce enough photons for the reioniza- 
tion while burning H and thus forming metals (Li and higher). 
These metals are released in supernovae explosions at the end 
of the stars short liv es (~ 1 Myr), where after they clump to- 
gether to form dust, iNozawa et a ll d2003h . Knowing the pro- 
duction rate of ionizing photons to be ^-/rit, (Fig. O, we can 
calculate the total photon energy released from the Pop III stars. 
This can be done by supposing that each photon has an effec- 
tive energy of E 7 — c y J v dv hvB v (T,)/ J y dvB v (T t ), where 
hvi on = 13.6 eV and B v {T t ) is the spectrum of a star with tem- 
perature T* . The energy of the non-ionizing photons is included 
through c y = u tot I u v>Vimi (u is the energy from the star). A Pop 
III star has T t ~ 80 000 K JShiova et alJl2002l page 9) which 
gives E y =s 36 eV. Note that for other reasonable star temper- 
atures, Ey does not vary significantly, Ey\ f>0xl0 3 K ~ 36 eV and 
40 eV. Hence, the total Pop III photon energy 
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production is Ey^/rij, per baryon per unit z- For each con- 
sumed nucleon, we assume that a nuclear energy of E r = 7 
MeV is released as radiation, which means that the nucleon 
consumption rate is j-'-jr/ni, nucleons per baryon per unit z- If 
fd is the fraction of the consumed baryon mass that becomes 
interstellar dust, (some of the metal atoms will remain with the 
core after the SN explosion, some will stay in the close vicinity 
of the SN and some will never clump together to form dust) the 
co-moving dust production rate will be 



Ey driy 

J+ = fd—®-b—rl n b- 

E r dz 



(4) 



A dust grain will eventually be destroyed, e.g. by col- 
lision, by supernova S hockwaves or by cosmic rays, see 



iDraine & SalpeteJ \ 19791) for further discussion. If a dust grain 
has a lifetime of At we can write the dust destruction rate as 



J = 



£^,o(z) 



Qrf,o(z) dt ^ 

At dz~ Ar// ( Ai/ 2 (l +z) 5 /2' 



(5) 



where £2 m is the relative matter content today, because the uni- 
verse is matter dominated for 5 < z < 15. 

Solving Eq.|3]gives the dust density evolution 



&d,o(z) 



(6) 



where z, = 20 is the beginning of the dust formation (see Fig.[Q 
and 



[2 (1 +ZT 3 ' 2 
Y(z) = exp\- —--T 

3 Q, l 2 AtH n 



(7) 



We note that the source term J + is modulated by the destruction 
term yf^-- The dust density is plotted in Fig. [5] where we note 
a strong dependen cy on the dust lifetime. In local dust At ~ 
100 Mvr. lDraine & Salpeteij i 19791) . Howe ver, the uncertainty 
is rather large, according to Draine ( 1990), Ar = 30 Myr - 10 
Gyr, depending on the environment. Note, however, that the 
density at the reionization red-shifts is much lower than in the 
interstellar medium in the Milky Way which implies a rather 
long dust life-time. 
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Fig. 3. The co-moving relative dust density evolution Qd,o = 
Pdusi/Pc, for fd = 1. The minima at z — 6 and 9 for At < 0.1 
Gyr is due to the fact that Az = 1 is not a constant time interval. 

3. Results and Discussion 

3.1. Metallicity 

If we suppose that most of the metals were ejected as dust (not 
as gas) the metallicity comes from the dust grains. The metal- 
licity is directly obtained through the produced dust. By letting 
At — > oo (At = 10 Gyr is good enough) we find the metallicity: 



Z 



Q rf ,u(Af -> oo) 

0.02 ■ n b 



1147 ■ Q dfi (At -> oo) 



(8) 



or in absolute terms Z « 22.90^0- At z — 5 we have Q^o = 
2.3 • l0- 5 f d , which gives Z « 5.2 • I0~ 4 f d = O.O26/ rf Z . 

There are not much metallicity data available for z > 5. 
Metal poor stars in our galaxy are one point of reference, ab- 
sorption lines in the Lya spectrum from quasars are another 
one. The lowest metallicities found in stars in the Milky Way 
are Z/ Z n ~ 0.01,|D epagne et al. (2002). The Lya forest sug- 
gests JSongaila & Cowidl2002l figure 13) that Z/Z G ~ 0.003 
for z ~ 4.5 assuming that [Fe/Hl « log(Z/Z Q ) as suggested 
by IVandenBer g et all EoOO. page 432). This indicates that 
fd ~ 0.1. However, th is might be lower than the actual value, 
cf. JPettinietalll997t figure 4). 

In heavy stars, virtually all the helium is consumed, pro- 
ducing metals. For simplicity (and lack of data), we assume 
that all the ejected metals clump to form dust, fy * f eject- This 
means that fd will almost entirely dep end on the dust ejec- 
tion rate in the supernova explosion. In llwamoto et alJ £[998) 
a detected hypernova of mass M ~ 14M Q seems to have 
f eA ert <; 0.7. Furthermore , according to a dust production model 
bv iNozawa et al.1 J2003I) . f d w 0.2 - 0.3. At the same time, 
some of the stars will become black holes, not ejecting any 
metals. lHeger & Wooslevl (12002), decreasing fy. Currently this 
decrease is largely unknown. 

In summary, the mass fraction of the produced metals in the 
Pop III stars, having become interstellar dust, should be around 
fd ~ 0.1 - 0.3. In the following we use the more conserva- 
tive f d — 0.1, in agreement with the Lya forest measurements, 
unless otherwise stated. 
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Fig. 4. Opacity r with dust evolution taken into account. 
3.2. Dust Opacity 

With our model for the dust density evolution, we want to cal- 
culate the opacity of the dust, as seen by the CMB. This will 
tell us how much the CMB spectrum is altered by the passage 
through the dust. 

The dust opacity is given by 



r dt 

= c I dz—o- v n d (z) 

Qoc 3pc [ d {— 



«*o(z)(l + z) 



(9) 



(10) 



where v (v e ) is the observed (emitted) frequency and v = 
v e /(l + z). The dust number density is tidiz) — (1 + z) 3 x 
p c Cld,o(z)/m g where m g = ^j-p g is the grain mass. We see 
(from t oc Q.d,o) that t is proportional to the parameter fd. 

The resulting opacity can be seen in Fig. |4] We note that 
the opacity is small, t « 1. The smooth knee is due to the 
change of j3 at the redshifted v r , see Fig. 12 but this is not in 
the spectral range of the CMB. The differential opacity dr/dz 
is plotted in Fig.[5]for A = 1 mm. We see that with a short dust 
lifetime, the dust density falls off almost immediately (in terms 
of z). However, for longer lifetimes, the early dust could still 
play a certain role for z < 3. This could eventually contribute 
to dimming of distant objects. We also note the impact of the 
expansion of the universe in decreasing the dust density and 
thus the opacity. This is why the increase in Fig.^ at z ~ 5, is 
not apparent in the opacity, Fig.|4] The submillimetre effective 
dust opacity follows a v 2 emissivity law. 

3.3. Dust Temperature 

In order to deduce the equilibrium temperature of the dust, we 
write the balance between the absorbed CMB, the absorbed 
starlight and the emitted IR light from the dust: 



Pd - P* + Pcmb- 

The powers P d and Pcmb can be written as 

P x =4tt dv e cT v B v XTx), 
Jo 



(11) 



(12) 
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Fig. 5. The differential opacity dr/dz at A — 1 mm for different ^ c j f - ■ ■ u t a, 

& v J ' b lg. 6. Energy density of ionizing photons compared to uqmb - 

Acts Ti MB /c, where as is Stefan-Boltzmann's constant. 



where B v is a Planck blackbody spectrum and X = {CMB, d). 
In the wave-length range considered, the spectral index B — 2. 
Supposing that B is constant, eq.^2 can be solved for the dust 
temperature analytically in the submm range: 



r 4 +J S 



T 4+P T 4+fl 
*eff CMB> 



where the effective temperature is defined by 



*eff 



(13) 



(14) 



and C B = £° dxx 3 ^ /(e x ~ 1) = + 3)! Zti k~ (A+fi \ such that 
C () * 6.494, d ~ 24.89 and C 2 ~ 122.1. 

However, in our calculations we use the exact eq. ^2 an d 
1121 while eq.^^can be used as a cross-check. 

The absorbed power density, P t from the radiation of Pop 
III stars peaks in the UV-region and can be approximated by 



(15) 



where cruv is the dust-photon cross section in the UV region 
and the energy density is 



w*(z) = fe 



X z , dn y 



l+z 
1+z' 



(16) 



where f esc is the escape fraction of photons from the star halos. 
We neglect the loss of photons due to the reionization itself. 
Eyjji is the effective energy of t he photon em itted at z' and 
then redshifted to z. According to ICen. Rj d2003l) . f esc = 0.3 
gives an electron opacity r e * 0.13 which is within one stan- 
dard deviation of the results by WMAP Hereafter, we adopt 
this value of f esc . 

The energy density of the ionizing photons are compared to 
the CMB in Fig. |6] The star energy density is much less than 
the CMB energy density at this epoch, and the curve resembles 
the accumulated photons in Fig. [2 Hence, the dust temperature 
closely follows the CMB temperature, see Fig.Qand Eg. 1121 
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Fig. 7. The dust temperature is plotted against the CMB tem- 
perature with the relative quatity (Tj - T C mb)/Tcmb- 

3.4. Observed Intensity 

Now we proceed to compute the average intensity (monopole 
term) of the submm and microwave background which is made 
of the CMB and early dust emission. The simple radiative 
transfer of the CMB through the uniform dust screen yields 
the following observed intensity: 



B v (Tcmb) + 



f 

Jo 



(17) 



(1 +zf 

From Fig.0]and0 we see that the opacity is small, (t «: 1) 
and the dust temperature is only slightly higher than the CMB 
temperature (7^ > Tcmb)- This gives the following formula for 
the excess intensity relative to the unperturbed CMB: 



Az'y = i v 
= Tcmb 



B v (Tcmb) 
dB v 
If 



T=Tcmb 



f 

Jo 



Td(z) - Tcmb (z) 
Tcmb(z) 



dr P 



(18) 



where Tcmb is the CMB temperature today. The integrant is 
plotted in Fig. Q We note that a new component is added to 
the primary CMB spectrum. Eq.^]tells us that it has a specific 
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Fig. 8. Comparison of the modeled intensity for the early dust 
emission in excess of the CMB with the observed FIRAS spec- 
trum (dash ed red curve) of the cosmic far IR background as 
detailed by Lagache et al1 dl999h . 



spectrum which is the product of a 2.725 K blackbody tem- 
perature fluctuation spectrum (like primary anisotropics) and 
a v 2 power law (from d.T e ). This effect is mostly visible in the 
submm range and has a minor contribution in the radio domain. 

In Fig.|8j the excess intensity is plotted alon g with the extra - 
galactic background measured by FIRAS, Pug et etalJ (fl996): 
iFixsen et alJ Jl998h : lLa gache et alJ lll999h . Depending on the 
dust destruction rate (parametrized by the dust lifetime Af), the 
computed early dust background can be an important part of 
the observed background from 400 /vm up to the mm wave- 
length. The exact position of A r will only slightly displace the 
spectrum, leaving the magnitude unchanged. Most of the far IR 
background can now be expl ained by a population of z — to 
z — 3 luminous IR galaxies, Gisp ert et all (2000). A fraction 
of the submillimetre part of this background could arise from 
larger redshift dust emission as suggested by Fig. [8] 

In order to check our results, we calculate the co-moving lu- 
minosity dens ity of the d ust in the submm region and compare 
it with ( Gispert et al. 2000, figure 4). We find them compatible. 



3.5. Discussion 

Just like the Thomson scattering during reionization, early 
dust will also tend to erase the primordial anisotropics in the 
CMB. However, due to the much smaller dust opacity (compare 
Td(lmm) < 10" 3 and r e = 0.17), this effect will be negligible. 

The early dust will also introduce a new type of secondary 
anisotropics with a typical size of a dark matter filament. Here, 
we only estimate the order of magnitude of this effect. If the co- 
moving size of the dark matter filament is L, the angular size is 
3 • (L/5 Mpc) arcminutes at z = 10 which corresponds to mul- 
tipole number t ~ 4000 ■ (L/5 Mpc). Fortunately, this region in 
/"-space does not contain any primordial fluctuations because 
of the Silk damping. However, there are ot her foregrounds in 
the same region, see Agha nim et alJ J2000t) . If we suppose a 
contrast of 10% in the dust intensity between dark matter fila- 
ments and the void, we obtain values of AT /T as 3 x 10~ 7 (for 



(2003) and FIRAS-II type of in- 
d200 



A — 1 mm, fd = 0.1 and Af = 1 Gyr). These anisotropics, 
pending more accurate calculations, clearly are in the range 
of expected arcminute secondary anisotropics from other ef- 
fects. They could b e detected by P lanck HFI (High Frequency 
Instrument), Lam arre 
strument. IFixsen & M ather (2002). 

The results of these calculations depend only very weakly 
on the precise dust model assumptions. We have also tried a 
different (but similar) shape of the ionizing photon production, 
Fig-ID an d found that the results do not vary significantly. 

Very little is known about the universe during the reion- 
ization epoch. Nevertheless, there are several parameters that 
could be calculated more accurately. 

The two most important parameters in the present model 
are the dust lifetime, Af and the mass fraction of the produced 
metals that are ejected as interstellar dust, fy. The dust life- 
time could be determined more precisely by making 3D sim- 
ulations of the dust production in combination with structure 
formation. The simulations would also give the inhomogeneous 
dust density evolution. The result would be a better estimate of 
the aforementioned secondary anisotropics caused by the vari- 
ations in the dust opacity. A more refined dust grain model, 
using e.g. a distribution of grain sizes would also be more real- 
istic. If the dust is long-lived, it could also have a certain impact 
on measurements in the optical and UV region. Finally, we note 
that most of the results are proportional to the dust density and 
thus to fj. To evaluate fy more precisely, we need a better un- 
derstanding of the typical properties of the first generation of 
stars, see section l3~T1 which is currently much debated. 

4. Conclusions 

We have shown that the radiation from early dust, produced and 
heated by Pop III stars, contributes to the extragalactic submil- 
limetre background within the limits set by FIRAS. It may not 
be detected by the present generation of instruments but fu- 
ture experiments such as Planck HFI and FIRAS-II should be 
able to measure it, by using its specific predicted spectral signa- 
ture. This high-redshift dust, contemporary to the reionization, 
should show up as small-scale anisotropics when observed by 
sensitive submillimetre instruments. These anisotropics are in 
the same range as other small-scale anisotropy effects. 
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